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^ ' We present Yang-Baxter maps associated to elliptic curves. They are related to discrete ver- 

sions of the Krichever-Novikov and the Landau-Lifshits equations. A lifting of scalar integrable 
quad-graph equations to two-field equations is also shown. 
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1 Introduction 

'<\ 

The importance of the Yang-Baxter (YB) equation to a variety of branches in physics and mathe- 
matics is well known. Its solutions are intimately related to exactly solvable statistical mechanical 
models, link polynomials in knot theory, quantum and classical integrable models, conformal field 
theories, representations of groups and algebras, quantum groups and many others. More interest- 
ingly, YB equation provides various connections among the aforementioned disciplines. 

Historically, YB equation has its roots in the theory of exactly solvable models in statistical 
mechanics \35\ [8] and the quantum inverse scattering method [31] . For an extensive account of early 
00 I work on the YB equation see jl4] . In its original form the quantum YB equation is the relation 

■ i2(2,3) ^(1,3) ^(1,2) ^ ^(1,2) ^(1,3) ^(2,3) (1) 

' in End^V^^), for a fc-linear operator R : V^"^ i-^ y®2^ where F is a vector space over a field k. 

Here, R^^'^^ is meant as R acting on the first and third factors of the tensor product V^^ and 
Q ! as identity on the second, and similarly for R^^''^^ and R^'^'^\ Drinfel'd suggested to study the 

simplest possible solutions of the YB equation by replacing (fc-Vect, (8") with (Set, x), where the 
YB equation is regarded as an equality of maps in for a finite set X. As it was pointed out 
in [TT], this setting provides potentially new interesting solutions of the original YB equation by 
I considering the free module generated by the set X. Various interesting examples of YB maps such 

as those arising from geometric crystalls have revealed a richer structure of the underlying 
set; X is an algebraic variety and R is a birational isomorphism. As in [33] we refer to solutions of 
the YB equation in (Set, x) simply as YB maps. One has to have in mind though that YB maps 
through a natural coupling can be regarded as equations on the edges of graphs. 

One of the most distinguished properties of integrable partial differential equations is their 
invariance under Darboux and Backlund transformations [21], [29]. The nonlinear superposition 
formulae of the solutions generated by the Backlund-Darboux transformations provide natural 
discrete versions of the continuous equations. In turn, the key transformation properties of the 
discrete equations are intimately related to the YB property for maps or its proper generalization 
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in higher dimensions namely the functional tetrahedron equation. A prime example of the latter 
are the star-triangle transformations in electric networks and the Ising model and their connection 
with discrete equations in the three-dimensional lattice associated with the Kadomtsev-Petviashvili 
(KP) hierarchy and its modifications [TSj. 

Recent interest on solutions of the YB equation for maps has appeared in the literature. This 
is mainly due to the development of the dynamical theory of YB maps [33], and classification 
results of YB maps for X = CP^ , in connection to analogous results for two-dimensional integrable 
discrete equations on the square lattice [5l[6]. The latter correspondence was further investigated 
in |26t l?7] by exploiting the local groups of symmetry transformations of the discrete equations. It 
was demonstrated that an integrable quadrilateral equation with a sufficient r-parameter symmetry 
group gives rise to a YB map. However, there exist discrete equations which do not admit any local 
symmetry group. Such an equation is a discrete version of the Krichever-Novikov (KN) equation 
[THl [IS], introduced by V. E. Adler in [3]. 

A characteristic feature of the discrete KN equation is that the lattice parameters lay on an ellip- 
tic curve. Remarkably, from the very early times exactly solvable two-dimensional models appeared 
in statistical mechanics it was observed that the solution of many of these models ultimately leads 
to the introduction of elliptic functions, such as the eight-vertex model on the square lattice [7J. 
Thus, it would be interesting to investigate whether there exist also solutions of the YB equation 
for maps related to elliptic curves. Already this problem was addressed in [21] where a theoretical 
framework was introduced for deriving YB maps from factorization of matrix polynomials and 
0-f unctions. 

The main aim of the present work is to exhibit YB maps with parameters living on elliptic curves 
and which are associated to integrable partial differential equations (PDE) . In Section [2] we present 
background material on the YB maps. In Section [3l we present key transformation properties 
of integrable lattice equations, encoded into braid type equations, and give a brief account on 
symmetry aspects of discrete integrable equations and their usage in deriving YB maps. Finally, 
we present a way for deriving a YB map from integrable lattice equations of certain type, without 
using a local symmetry group. The latter method is applied to generic integrable lattice equations 
associated to elliptic curves such as discrete KN equation and discrete Landau-Lifshits equation 
and the results are presented in Sections H] and [Sj respectively. The paper concludes in Section [6] 
with various comments and perspectives. 

2 Definitions and notation 
Braid and YB maps 

In the following we use the notation and terminology introduced in [33j (for a recent review see 
j34| ) . Let X be an algebraic variety, and iZiXxX^XxXa birational isomorphism. Let 
R{i,j) : X" ^ X" denote the map acting as R on the components of the n-fold Cartesian 

product X X X X • • • X X and as the identity on all others. More explicitly, for y G X let us write 





Then, for n > 2 and 1 < 



i, j ^ n, i ^ j the map R^'^'^^ is given by 




, Xj^i, QyXi , Xj ) , Xj-\.\ , . . . Xfi ) , 



Xi—\i f{Xi , Xj ) , Xi-i-1 , . . . Xji) ) i ^ j • 



i < j, 



(3) 
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In particular, for n = 2 we have = R and R^'^^^\x,y) = {g{y,x)J{y,x)). The latter map IS 

R conjugated by the permutation map a, defined by a{x,y) = {y,x), i.e. 



R^'''^=aRa. 

Definition 2.1. (i) The map R is called a YB map if R satisfies the YB equation 

^(2,3) ^(1,3) ^(1,2) ^ ^(1,2) ^(1,3) ^(2,3) ^ 

regarded as an equality of maps of into itself. 

(a) R is called reversible, or unitary, if it satisfies the condition 

R^^'^^ R = ld^2. 



(4) 



(5) 



(6) 



(Hi) R is called non- degenerate if the maps from X into itself defined by s ^ f{s,y) and 
t g{x, t) are bijective rational maps for any fixed y and x, respectively. 

A schematic representation of the YB equation is given by the two decompositions of an ele- 
mentary 3-cube as depicted in figure [H The composition of maps in the LHS and RHS of the YB 
equation ^ are given by 



(a) : {xi,X2,X3) > {X^,X2,X3) > {Xi,X2,Xr^) y (Xi,X2,X3) , 

/j(2,3) RC'--'^') _R(^'2) 

(b) : {xi,X2,X3) — > {xi,X2,xl) — > {xl,xl,xl*) — > {xl*,xl*,xl*) , 



(7) 



respectively. The YB equation guarantees that the images of (xi,X2,X3) G X^ under the two 
composition of maps in ^ are identical, thus the two parts of the 3-cube can be glued together. 





(a) ■ (b) 

Figure 1: A cubic representation of the YB relation 
Remark 2.2. Defining B = aR the YB equation ([5]) for R translates to the braid equation 

^(2,3) ^(1,2) ^(2,3) ^ _g(l,2) ^(2,3) ^(1,2) ^ 

for B. The unitarity condition ([6]), in view of takes the form 

B^ = M^2 , 



(8) 



(9) 



and B is said to be an involution. 



3 



Lax matrices for YB maps 

Instead of a single map one may consider a whole family of YB maps parametrized by two continuous 
parameters ai, a2 G C, where C is an algebraic set in C^. The YB relation then takes the parameter- 
dependent form 

p(2,3) p(l,3) p(l,2) _ p(l,2) „(1,3) p(2,3) . ^ 

-ft/ ^ It, S -ft/ N -ft/ \ -ft/ N -ft/ N , IJ-UI 

(02, as) (01,03) (01,02) (01,02} (01, as) (02,03) ^ ' 

and the unitarity (reversibility) condition becomes 

,= Idj'z . (11) 
(02,01) (01,02) ^ ^ ' 

In the following we drop the dependence on the parameters and write just R for a two-parameter 
YB map, since we always consider maps of this type. 

Let L{x; a, A) € Mat(r, X) be a two-parameter family of r x r matrices depending on x € X and 
polynomially/rationally on the coordinates of a, A S C. The following notion of Lax matrix for a 
YB map was introduced in [30]. 

Definition 2.3. (i) L{x;a,X) is called a Lax matrix of the YB map R, if the relation R(xi,X2) = 
{xi,x'2) implies that 

L{x2] 02, A) L{xi;ai,\) = L{x'i;ai, A) L{x'2; 02, A) , (12) 

for all A € C L{x] a, A) is called a strong Lax matrix of R, if the converse also holds, 
(ii) L(x;a,X) satisfies the n-factorization property if the identity 

L{x'n;an,\) ■ ■ ■ L{x'2]a2,X) L{x^;ai,\) = L(x„;q„,A) • • • L(x2; 02, A) L(2;i; ai, A) , (13) 

over C, implies that x/ = Xi, i = 1, . . . , n. 

Remark 2.4. The 2-factorization property of L corresponds to the unitarity property of R, while 
the 3-factorization property to the YB property. Indeed, the composition of maps 

(X1,X2) {x'l,x'2) {x'lx'i), (14) 

is represented by the matrix factorization 

L{x2; OL2, A) L(xi; Qi, A) = Lix'^, qi. A) L(x2; 02, A) = L{x2\ 02, A) L{x'[\ ai. A) . (15) 

Evidently, the unitary property of R is equivalent to the 2-factorization property of L. On the other 
hand, the cubic representation of the YB relation (see figure [T]) suggests to consider the product 

L(x3;a3,A)L(x2;a2, A)L(xi;ai, A) . 

It can be factorized in two different ways, according to the composition of maps in ([7|), i.e. 

(a) L(x3;a3,A)L(x2;a2, A)L(xi;ai, A) = L(x3; 03, A)L(x'i; oi, A)L(x2; 02, A) 

= L{x'[\ ai, A)-L(x3; 03, A)L(x2; a2, A) 
= L[x'[\ ai, A)-L(x2; 02, A)-L(x3; 03, A) , 

(b) -L(x3;a3,A)-L(x2;a2,A)-L(xi;ai, A) = L(x2; 02, A)L(x3; 03, A)-L(xi; ai, A) 

= L(x2; a2, y)L{x\\ ai, A)L(x3*; 03, A) 
= L(xi*; ai, A)-L(x2*; a2, A)L(x3*; 03, A) . 

where we have used the fact that L is a Lax matrix for R in each face of the cube and the associativity 
of matrix multiplication. Thus, the YB property of R is equivalent to the 3-factorization property 
of L. 
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Example 2.5. Consider the map R : C 

R{x,y) 



X CP 



X CP^ defined by 



ai+xy a2 + xy 

y ) X 

a2+xy ai+xy 

introduced in [27J. The above map admits the Lax matrix 



(16) 



L{x; a, A) 



X 
X 



2 \x 



a 



It is straightforward to check that the discrete zero curvature equation (jl2p . imphes the map (jlOp . 
Conversely, according to [30] a hint for considering the matrix L of the above form, is based on 
the YB map itself. Indeed, one notices that the second component of the map can be written as a 
linear fractional transformation induced by the linear transformation with matrix L[x; ai, X)\\=a2 
on [y 1]^. For the above Lax matrix the n-factorization property can be proved as follows. The 
null space of the linear transformation in the LHS of (|13|) . for A = ai, is spanned by the vector 
[— ai x'^ . Similarly, [—ai xi\" spans the null space of the RHS linear transformation. Because 
of the identity (jl3p . we conclude that xi = x'l^ and the rightmost matrices cancel out. Therefore, 
by induction, L satisfies the n-factorization property. The YB map ()16|) is simply related to the 
Fjjj map obtained in the recent classification [6j for the case X = CP^. 



3 YB maps and integrable quadrilateral equations 

In this section we present first braid transformation properties of integrable discrete equations 
defined on elementary squares. Next, we briefly summarize a method for obtaining YB maps from 
integrable discrete equations on quad-graphs, which is based on the existence of a local group of 
symmetry transformations of the equations. Finally, we present another method to the same end 
which does not prerequisite the existence of a local group of symmetry transformations. 

Main properties of integrable discrete equations 

We consider discrete equations on quad-graphs given by an algebraic equation 

Q(/A,/B,/r,/zi;a,/3) =0, (17) 

relating the values of a function f -.1? assigned on the four vertices of an elementary plaquette. 
It is assumed that (i) opposite edges on the plaquette carry the same lattice parameter a, j3 and 
(ii) equation (fT7|) it can be solved uniquely for each /j, say /r, i.e. 

fr = HfAjBjA;a,(3). (18) 

In order to make contact with the special properties of the integrable discrete equations we interpret 
equation (I17p as a map jB : — > defined by 

B{fAjAjB) = ifAjrjB), (19) 

where fr is given by (fTSl) (see fig. ([2]) ). 
Let us now define Bj : X" ^ X" by 

Bj = ldxx ■■■ X B X ■■■ xldx, (20) 

where B acts on the j — 1, j and the j + 1 factors of X" with parameters (aj_i,aj). The key 
properties of maps associated to integrable discrete equations on quad-graphs are the relations 

B] = Idxn , {Bj Bj+if = Idxn , Bj B^ = BiBj, |i - j 1 > 1 , (21) 
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Figure 2: An elementary quadrilateral for a quad-graph equation. The arrow indicates the flip 
{fA,fA,fB) ^ {fA,fr,fB)- 

see [1]. The first one means that each transformation Bj is an involution. The second one, in view 
of the first, yields the following braid-type relation 

Bj+iBjBj+i=BjBj+iBj. (22) 



A' A' 




r r 



Figure 3: A cubic representation of the braid relation (|22p for quad-graph equations. 

The braid relation ()22l) guarantees that the three-dimensional consistency property [23j . [9], 
which nowadays is synonymous with the integrability of a quad-equation holds (for a recent account 
on the subject we refer to the monograph |10j). 

Prom integrable discrete equations to YB maps via symmetry groups 

Local symmetry groups of transformations of integrable discrete equations provide a natural way for 
obtaining YB maps from them. The main observation is that the variables of certain YB maps can 
be chosen as invariants of the symmetry group admitted by the corresponding lattice equation. The 
symmetry approach was exploited in [26], where it was also shown that all classified quatrirational 
YB maps, for X = CP\ found in [6], can be constructed from integrable quadrilateral equations. 

Definition 3.1. Let G be a one-parameter group of transformations on X'^, of the form 

G: {x,y,z)^iX{x;e),Yiy;e),Ziz;e)) , egC. (23) 

and B a map o/X'^ into itself. G is said to be a local (Lie-point) group of symmetry transformations 
of the map B ifGoB = BoG, for every e € C. 

Let X = CP"*^ and consider the following map 

B{x,y,z) = {x,y + - — -,z), (24) 
X — z 

which is associated with the discrete KdV equation. The corresponding map Bj defined by ([^U]) sat- 
isfies the braid type relation ()22p . Moreover, the map (I24p commutes with the group of translations 
given by 

Gi: {x,y,z) ^ {x + e ,y + e ,z + e) . (25) 
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Thus, Gi is a Lie-point symmetry of the map (j24p . The action of Gi on is regular with one- 
dimensional orbits, thus local coordinates on the set of orbits of Gi are provided by the complete 
set of functionally independent invariants for the group action: 



u = y — X , V = z — y . 
Projecting the map (|24p to the set of orbits of Gi we obtain the map 



B{u,v)={u-^,v + ^] , (26) 



U + V U + V 

which satisfies the parameter braid relation ([H]). Thus, the map 

^, ^ ( a — b a — b\ 
R{u,v)=aB{u,v)= (v + ^—,u-^—) , 27 

\ U + V U + V J 

is a YB map, known as the Adler map [2]. The most general local group of symmetry transforma- 
tions of the map is G = 5*0(1, 1), generated by Gi and the one-parameter subgroups G2, G3 
given by the group actions 

G2:ix,y,z) ^ {x-e2,y + e2,z-e2), 
Gs:{x,y,z) (x 6"^=^ , y e''^ ,2:e"^^). 

By using similar arguments one may consider the set of orbits of the subgroups G2, or G3, to obtain 
other YB maps from the map (j24p . More precisely, we have the following 



Proposition 3.2. Let X = and a map B :'K^ ^ li^ satisfying the braid-type relation ^22\}. If B 
admits a local group G of symmetry transformations which acts regularly on with n- dimensional 
orbits, then the projection of the map B to the set of orbits of G satisfies the braid relation. 

Proof. The assumptions for the action of G on X^ guarantee the existence of a 2n-dimensional 
quotient manifold denoted by X^/G, i.e. the set of all orbits of G. Local coordinates {u,v) G X^ 
can be chosen by a complete set of functionally independent invariants for the group action, see 
e.g. Theorem 3.18 in [25]. 

Let us denote by S : X^ 1— > X^ the projection of the map B on X'^/ G. The braid property of the 
map B is inherited by the braid type relation ()22p which satisfies the map B. This can be easily 
deduced from the cubic representation of the relation (|22p (Figure E]) . It should by noted that the 
invariants of the group action (YB or braid variables) can be naturally assigned to the edges of the 
elementary squares instead of the vertices where the variables of the original map B are assigned 
to. □ 



Two— field integrable discrete equations as YB maps 

The existence of a symmetry group of the map B provides us a way to obtain a YB map by using 
as YB variables the invariants of the group action which are naturally attached on the edges of 
the squares. As it was shown recently [28], [32], generic integrable quad-graph equations, such as 
the KN discrete equation, do not admit any local symmetry group of transformations. Thus the 
question arises whether such equations are related to YB maps, as well. This question is answered 
in the affirmative in section [H Proposition 13.31 below shows how to cast two-field quad-graph 
equations of a certain type into YB map form. Moreover, it motivates a way of lifting an integrable 
scalar quad-graph equation to a two-field one and consequently to recast the equation into a YB 
map. 
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Specifically, we consider lattice equations where at each vertex there is a two-field (u, v) G X'^ 
and the defining relations on the quadrilateral are (see figure [2]) 



{ur, vr) = {Fi{uA,UB,VA;a,b) , F2{vA,UB,VA;a,b)) 



(28) 



where Fi,F2 take values in X. This scheme of two-field quad-graph equations, although not the 
generic one since it does not involve all eight values of the fields, arises in the superposition formulae 
of Backlund transformations for two-field integrable PDEs e.g. the nonlinear Schrodinger system 
[TB] , PP . The aim is to recast discrete equations of the form ([^5]) into a YB map form. To this end 




Figure 4: The quadrilateral with "thickened" edges as a parallelepiped and the assignment of the 
YB variables for the map ()33p . The variables xi and qi are identified, respectively y2 and p2, 
(6-point scheme). 



we group the fields appearing in the RHS of equation (j28p as follows 

{xi,x2) = {ub,va) , (yi,y2) = {ua,va) ■ 
A pictorial representation of the assignment of the YB variables for the map 

R : {{xi,x2), (yi,y2)) ((pi,P2), {qi,q2)) , 
(equation ()33p below) is shown in Figure U] where 

{pi,P2) = {ur,VA) , (91,92) = {uB,vr) ■ 

Proposition 3.3. Let the map S : X^ x X^ x X^ ^ x x X^ defined by 

B{{ua,va),{ua,va),{ub,vb)) = {{uA,VA),iur,vr),{uB,VB)) , 

where {ur,vr) is given by /i28]). satisfies the braid relation /i22\) . Then the map R : X^xX^ 
defined by 

R{{xi,X2), (^1,2/2)) = {{Fi{yi,xi,y2;a,(3) , 2/2) ,{xi , F2{x2,xi,y2;a, (3))) , 
satisfies the YB relation. 



(29) 



(30) 



(31) 



(32) 



[2xX2 



(33) 



Proof. By straightforward calculations, we derive first the relations for the functions Fi, F2 such 
that B satisfies the braid type relation. The values {ur'-,vr>) and {ua>-,va>) are found in two 
different ways, according to the left and right hand side of the braid relation and are given by 

UK = Fi{ur,UA,VB;b,a) ua' = Fi{uA,UE,VK;c,a) = Fi{uK,UA',VB;c,b) 

VK = F2{vr,UA,VB;b,a) va' = F2{vA,UE,VK;c,a) vp' = F2{vk,ua' ,VB;c,b) 

UK' = Fi{uA,UE,vr;c,b) up' = Fi{ur,UK',VB;c,a) ua' = Fi{uK',UE,vr';b,a) 
VK' = F2ivA,UE,vr;c,b) vp' = F2{vr,UK' ,VB]c,a) va' = F2{vK',UE,vr';b,a) 
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respectively. Thus, we have the following functional relations satisfied by Fi, F2: 

Fi(Fi(ur, UA, vb; b, a), Fi{ua, ue, F2{vr, ua,vb; b, a); c, a),VB;c, b) = 

Fi{ur,Fi{uA,UE,vr;c,b),VB;c,a) , (34) 

F2{F2{vr, ua,vb; b, a), Fi{ua,ue, F2{vr,UA,VB; b, a); c, a), f^; c, b) = 

F2{vr,Fi{uA,UE,vr;c,b),VB;c,a) , (35) 

Fi{uA,UE,F2{vr,UA,VB;b,a);c,a) = 
Fi{Fi{uA, UE, vp; c, b),UE,F2{vr,Fi{uA,UE, vr; c, b),VB,c, a),b, a) , (36) 

F2{vA,UE, F2{vr, ua,vb; b, a); c, a) = 
F2{F2{vA,UE, vr; c, b),UE, F2{vr, Fi{uA,UE,vr; c, b),VB;c, a); b, a) . (37) 

On the other hand, the YB relation for the map (I33p gives the following functional relations for 
Fi, F2 

Fi{Fi{zi,yi, Z2; f3,-/), Fi{yi,xi, F2{y2,yi, Z2; I3,j)] a,-f), Z2; a, P) = 

Fi{zi, Fi{yi,xi,y2;a, /3), Z2;a,-f) , (38) 

F2{F2{y2,yi,Z2;l3,-f),Fi{yi,xi,F2{y2,yi,Z2;P,"f);a,j),Z2;a,P) = 

F2{y2, Fi{yi, xi,y2; a, f3), Z2; a,-f) , (39) 

Fi{yi,xi,F2{y2,yi,Z2;P,-f);a,-f) = 
Fi{Fi{yi,xi,y2;a,(3),xi,F2{y2,Fi{yi,xi,y2;a,P),Z2;a,-f);l3,j) , (40) 

F2{x2,xi,F2{y2,yi,Z2;P,-/);a,-/) = 
F2{F2{x2,xi,y2;a, /3),xi, F2{y2, Fi{yi,xi,y2;a, f3), Z2;a,'y); f3,-f) , (41) 

and two additional equations which are trivially satisfied. 
Making the following substitutions 

Up zi UA^ yi Ue ^ xi a 1-^ J, 
Vb ^ Z2 Vr ^ y2 va ^ X2 6 1— > /3 

c I— > a 

in equations ([M|) - ([37|) . the latter become identical to ([38|) - ()^T]) . respectively. □ 
Remark 3.4. Consider the case Fi = F2 = F i.e. 

{ur , vr) = {F{uA,UB,VA;a,b) , F{vA,UB,VA;a,b)) . (42) 

If Ui = Vi, i = A,B,A, then from equation (fl2|) we have ur = vr and the map (f32]l essentially 
reduces to a single field map, namely 

B^{ua,ua,ub) = {uA,ur,UB) , (43) 

and B can be thought as a lift of This observation suggests to lift the discrete KN equation to 
a two-field quad-graph equation and then write it as a YB map. The lifting process can be applied 
to all scalar integrable quad-equations listed in [5]. 
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Example 3.5. The simplest equation of the classification in [5j is the discrete (potential) KdV 
equation, namely 

/n+l,m+l — fn,m H 7 7 ; (^4) 

Jn+l,m Jn,m+1 

where {n,m) S Z^. Its lift obtained by equation (I42p takes the explicit form 

a — b a — b 

Un+l,m+l = Un,m H ^ , Vn+l,m+l = 'yn,m H — , (45) 

and satisfies the braid-type relation (|22|) . The corresponding YB map obtained by using Proposition 
13.31 reads 

{Pl,P2) = {yi + — ,^2), {qi,q2) = [xi, X2 + — ). (46) 

V xi-y2 J V xi-y2j 

The YB map (|46|) was derived in [17J from matrix factorization and is symplectic with respect to 
a canonical structure. On the other hand, equations (I45p are the Euler-Lagrange equations for the 
discrete variational problem associated to the following Lagrangian density 

C = Un+l,m Vn,m " Un,m Vn,m+1 + (a - 6) In {Un+l,m - Vn,m+l) ■ (47) 

The problem of the Lagrangian formulation of the quad-graph equations classified in [5J has been 
addressed recently in |20j . 

Remark 3.6. The YB maps ()33p obtained by this method are non-quadrirational. The quadrira- 
tionality property of maps as was introduced in [6] is equivalent to the nondegeneracy property 
which is often imposed additionally on the YB maps. This can be seen immediately since the map 
(xi,X2) {Fi{yi,xi,y2;a, P) , 1/2) is independent of X2. 

4 Lifting discrete KN equation to a YB map 

The master scalar integrable quad equation listed in [5], in the sense that the rest integrable 
discrete equations can be derived from it by proper degenerations of the elliptic curve or limiting 
procedures, is discrete KN equation [3]. Using the identification (/^, /s, /r, /z\) = {x,y,w,z) on 
the quadrilateral (Figure E]) the latter equation reads 

aB — bA , s 

a (xy + wz) — b (xz + wy) TPrixw + yz — ab[l + wxyz)) = . 

1 — a^o^ 

This is the form introduced by Hietarinta in where the parameters a = (a, A) and b = (6, B) 
lay on Jacobi quartics given by 

£ = {{a,A)eC^ : A^ = a^ + ka^ + l} , 

and k is the modulus of £. The binary operation defined by 

'aB + bA {AB + kab){l + a%'^) + 2ab{a'^ + 6^) 



a 6 : 



1 - a262 ' (1 _ a252)2 



endows the set £ with an abelian group structure, in which e = (0, 1) is the identity element and 
the inverse of a point a = (a. A) is the point —a = (—a. A). In the following we use the notation 
X = {x,X) etc, for points in C^. 
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Proposition 4.1. The map R{x,y) = {p,q) defined by 

{p, P) = {F{y, X, Y; a, b) , Y) , {q, Q) = {x , F{X, x, Y; a, b)) , 

where 



(48) 



F{x,y,z;a,b) 



{I - a? lf){b z - ay) X + {a B - b A){y z - a i 



{aB- bA){abyz- + (1 - a'^h'^){az- by) ' 
is a unitary YB map, with Lax matrix given by 

L{x; a, A) = p{x; a)^^^"^ W{x; a, A) , 

where 



(49) 



W{x;a,\) 



-XX -X 



a xX + \- 



aA - AX 
1 - a2A2 

aA - 



a 1 + XxX 



Xx + X 



aA - AX 
1 - a2A2 

aA - ^A 
1 - a2A2 



1 - a2A2 
and the scalar function p is 

p{x;a) = {x^X^ + l) a'^ -x^ -X^ + 2AxX , 

a,X e £. 

Proof. First we prove that the matrix L{x; a, A) is a Lax matrix for the map (j48p by showing that 
the factorization problem 



L{y; b, A) L{x; a, A) = L{p; a, A) L{q; b, A) , 
is equivalent to equations (jiSl) . Taking into account that A G f we have 

ELoEl=oS^M^,y■,a,b) X^K' 



(50) 



(a2A2-l)(62A2-l) 



{W{y-b,X)W{x-a,X)^ 



Equating the different powers of A, A, the matrix equation (jSOp is equivalent to the following system 
of algebraic relations 

Sijkeix,y;a,b) _ Sijki{q,p;b, a) 



p{x; a) 1/2 6)1/2 „)i/2 ^(g. b)i/2 ' 

i,j = 1, 2, A; = 0, 1 and ^ = 0, . . . , 6. We calculate the terms 

Snooix,y;a,b) = abx{y - X) , 
Siwo{x, y; a, b) = a^b^Y {X -y), 
Si2Qiix,y;a,b) = ab{X - y) , 

5'i2io(a^i b) = abX y {bY — ax) + bx — aY + a B y — b AX , 
'S'2iio(^C5 y, t>) = Xy{bY — ax) + ab{bx — aY) + xY{aBx — bAy) . 



Sijke{x,y;a,b) _ Sijke{q,p;b,a) 



(51) 



From system (j5T| we get 



Si2oi{x,y;a,b) Si2oi{q,p;b,a) 
For (ijke) = (1100) and (ijki) = (1160) equations ([52]) lead to 



(52) 



q = x, 



P = Y, 



(53) 
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respectively and using them, equations (|52p for {ijk£) = (1210), (2110) lead to a linear system 
which is uniquely solved for {p,Q), yielding 



p = F{y,x,Y-a,b), Q = F{X, x,Y; a,b) 



(54) 



With the solution given by (j53p . ()54p by straightforward calculations we find that system (j5ip is 
satisfied. 

Next we prove that the Lax matrix L given by ()49p satisfies the n-factorization property. For 
A = a we obtain 

W{x;a,a) = a ^ [-X 1] . 

Thus the kernel of the linear transformation 

L{xn; an, ai) • • • L{x2; a-^, ai) L{xi;ai, ai) , 

is spanned by the vector [1 Xi]'^ which leads us to conclude that Xi = X[ in ([TB]) . Likewise, for 
A = e we have 



W{xn; an, e) • • • W{x2; a2, e) W{xi;ai, e) 



n «i 



Xi) 



.i=2 





1 




) 


Xn 


[-Xi 



In this case the kernel of the linear transformations in the LHS and RHS of (jl3p is spanned by 
the vectors [1 x[]'^ and [1 xi]'^, respectively. Thus, xi = x\ and consequently L(x\;a\,X) = 
L{xi; ai, A) for all X € £. Hence, the number of matrices in equation (fT3]) is reduced by one and by 
induction the n-factorization property is proved. Finally, by remark [27il the Proposition is true. □ 



5 A discrete Landau-Lifshits equation as a YB map 

In the literature, there exist several discrete versions of the Landau-Lifshits equation representing 
the nonlinear superposition formula for the solutions generated by the Backlund auto-transformation 
[22] , [I], [1] . Here, we use the one introduced in [1] and we present the end result, namely the 
corresponding YB map derived from the lattice equations by using Proposition 13.31 The map reads 
the form 

R{{x,X),{y,Y)) = {{Fi{y,x,Y;a,b),Y), {x, F2{X,x,Y;a,b)) , (55) 



where 



Fi{x,y,z;a,b) 



K{y,z)x- L{y,z) 
M{y,z)x + N{y,z) 



F2{x,y,z;a,b) 



K{y,z)x + L{y,z) 
-M{y,z)x + N{y,z) ' 



Kiy,z)-N{y,z) 
K{y,z)+N{y,z) 
L{y,z) 
M{y,z) 



2s2yz- {asi + S3) {y - z) - 2as2 - 4/3si , 
(y -I- z) (a h{asQ -|- 3^2) + 4/3si -|- 3qs2 + S4) /(« - ^) > 
S3t/z -I- {as2 + 2/3si)(y - -z) + 4/3s2 - o? s\ , 
s\yz^S2{y- z) - S3, 



Ba 



m—l 



+ Ab 



m—l 



The parameters a = {a, A), b = (b, B) lay on the WeierstraB elliptic curve 



(56) 
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where a, /3 are complex constants, the invariants of the curve. It should be noted that in contrast 
to the previous case the functions Fi , F2 are different reflecting the fact that the discrete equations 
constitute a genuine two-field system. The Lax matrix introduced in [ij, is also a Lax matrix for 
the map (j55p and is given by 



L{x;a,X) = p{x;a) ^^"^ W{x; a, X) . 



(57) 



Here, 



p{x; a) = 2 A{r'^ + x - X + a) , 
the matrix components of W{x; a, A) are given by 

{Wix;a,X))^^ = {A + A)r + {X-a)iX + a-X), 

{W{x;a,X))^^ = {A + A){X + x-X)-{X-a){X + 2a)r-2A{\-a), 

{W{x;a,X))^^ = A + A-{X-a)r, 

{W{x;a,X))^^ = -{A + A)r-{X-a){X + a + x), 



and 



— {xX + a{x- X) +a + 2a^) . 



The Lax matrix ()57p satisfies the n-factorization property for n = 2,3. Indeed, first we note that 
for A = Ao, where 



Xo = -2a-A\ Ao 
the matrix W takes the dyadic form 



-A - A'^{Xo- a) , A 



3a2 + a 
2A 



W{x;a,Xo) = (Aq - a) 



■{2A- A{a- x)) 
(a - x){a + X) 



a + X 2A-A{a + X) 



Thus, the kernel of the linear transformation with matrix L(a;;a,Ao) is spanned by the vector 
[A{a + X) — 2 A a + X] from which we conclude that X = X' . Next, for general A, the product 
of two Lax matrices takes the form 

1 ^ ^ 

{L{y- b, A) L{x; a, X)^^ = E E ^^M^^ ' 

pv I HKyj k=Q e=o 

where we have used that (A, A) G £. Equating the different powers of A, A, the matrix equation for 
the 2-factorization is equivalent to the system of algebraic relations 



Sijke{x',y';a,b) 



SijM{x,y;a,b) 



p{x'; a)i/2 b)i/2 a)i/2 p(^y. 5)1/2 • 

For (ijki) = (1221), (1211), (1120) we find that 

5*1221 = -{X + y), S1211 - 5ii2o = -x{X + y) . 

Hence, x = x' and consequently L satisfies the 2-factorization property. 
Likewise, the product of three Lax matrices reads 



(L(z; c, A) L(y; 5, A) L{x; a, X)^ 



Efc=o ELo Sjjkejx, y, z; a, b, c) X^ A^ 
p{x;ay/^p{y;by/^piz;cy/^ 
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For (ijke) = (1241), (1221), (1130) we find that 



5i24i= {X + y){Y + z), Si22i-Su30 = x{X + y){Y + z). 

From the corresponding ratios we have that x = x' , and the number of matrices is reduced by one. 
The 3-factorization property is reduced to the 2-factorization property, which is satisfied, and the 
map (|55|) is a unitary YB map. 

6 Conclusions 

Two famiUes of YB maps with parameters hving on elhptic curves are presented. Both of them 
are based on the combinatorics and the geometry of a certain type two-field quad-graph system 
(6-point scheme) that allows to cast them into YB map form. It is this scheme that suggested the 
lifting of scalar integrable quad-graph equations to two-field ones and subsequently the derivation 
of their YB form. 

We end by giving a rough account on YB maps in x arising from two-field integrable 
quad-graph equations. In [27j such YB maps were derived by exploiting the symmetry groups of 
the equations listed in [IJ. In the present work it is shown how all discrete equations listed in [I] 
are casted in YB map form (Proposition 13. 3j) . This list of YB maps is enhanced by "lifting" all 
integrable quad-graph equations listed in [5J , as it was demonstrated here for the generic equation 
of the class, namely the discrete KN equation, denoted by Q4 in [5] . Moreover, the list of YB maps 
is enriched by considering the symmetry groups of the lifted discrete equations. Thus, it turns out 
that even in this particular case (corresponding to the 6-point scheme) one has already quite an 
amount of YB maps in x and the problem of their classification becomes interesting in order 
to: i) find the representatives up to equivalence with respect to some group of transformations and 
ii) make the list exhaustive. 
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